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SUMMARY 

The  drag  of  artificial  earth  Bate llite 5  is  reconsidered  in  the  light  of 
recent  studies  of  gas-surface  interactions  and  atmospheric  composition. 
Between  heights  of  140  and  400  km  at  times  of  low  solar  activity,  or  heights 
of  140  and  600  km  at  times  of  high  solar  activity,  the  drag  coefficient  is 
almost  independent  of  height,  and  at  present,  because  of  the  lack  of  decisive 
evidence,  there  is  perhaps  not  sufficient  reason  to  abandon  the  value  of  2.2 
which  has  been  widely  used  in  recent  years  for  the  drag  coefficient.  It  oust 
be  recognized,  however,  that  this  value  is  subject  to  some  uncertainty  and  may 
be  too  low,  perhaps  by  as  much  as  10  per  cent.  At  heights  above  A 00  km  (leer 
solar  activity)  or  600  km  (high  solar  activity)  the  drag  coefficient  increases 
with  height,  since  both  the  degree  of  energy  transfer  and  the  moiawilar  speed 
ratio  decrease  as  the  molecular  weight  of  the  atmosphere  decreases. 
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1  IN'IRCDUOTICN 

Sinoe  a  satellite  normally  experiences  free-molecule  flew,  knowledge  of 
the  mechanism  of  molecular  reflexion  at  its  surface  is  necessary  if  the  aero¬ 
dynamic  drag  is  to  be  determined.  When  the  early  satellites  were  launched 

only  meagre  information  on  gas-surfaoe  interactions  was  available;  drag  was 
1 

estimated  v assuming  the  atmospheric  molecules  were  diffusely  re-emitted  with 
velocities  corresponding  to  the  surface  temperature  of  the  vehicle.  Even  today 
there  is  still  no  experimental  evidence  on  surface  phenomena  in  the  appropriate 
energy  range,  although  recent  theoretical  studies  have  indicated  the  important 
factors  involved.  In  the  light  of  these  investigations  it  is  new  possible  to 
indicate  how  the  drag  coefficients  may  differ  from  the  values  currently  in  use. 

2  jftEB-kCIECULE  FLCV? 

The  conventional  requirement  for  free-molecule  flow  is  that  the  mean 
free  path  in  the  free  stream,  Xq,  is  large  compared  with  a  typical  linear 
dimension  L  of  the  body.  The  ratio  X  is  known  as  the  Knuds en  number  K,  so 
that  free-molecule  flew/  occurs  if  K  »  1.  Far  a  body  moving  at  very  high  speech, 
however,  this  condition  is  not  sufficient  to  ensure  free-molecule  flow  if  the 
surface  temperature  is  low  and  appreciable  energy  accommodation  takes  place. 

If  the  number  density  in  the  free  stream  is  nQ  and  the  molecules  impinge 

on  the  body  with  a  speed  v^,  the  number  of  molecules  arriving  at  unit  projected 

area  in  unit  time  is  n  v,.  The  molecules  will  be  re-emitted  at  a  speed  v 

o  i  r 

which  generally  differs  from  v^.  The  number  density  of  re-emitted  molecules 

near  the  surface  will  then  be  of  the  order  n  v./v  •  Since  mean  free  path  is 

oiri* 

inversely  proportional  to  density,  the  mean  free  path  of  an  incident  molecule 
passing  through  the  re-emitted  molecules  near  the  surface  is  given  by 


If  free-molecule  flow  is  to  exist,  the  Knuds  en  number  based  on  X  oust  be  greater 
than  unity  and  the  true  requirement  is  X/L  »  1,  i.e. 


(D 


In  practice  the  largest  value  of  v^/vf  for  s  satellite  is  about  13,  but  this 
will  only  occur  when  the  molecules  are  re— omitted  with  a  speed  corresponding 
to  the  surface  temperature. 
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Values  of  mean  free  path  given  by^  the  U.S.  Standard  Atmosphere  (1962)  are 
sufficiently  accurate  for  order-of -magnitude  studies  and  are  shown  in  Fig.1  for 
heights  between  100  and  250  km*  At  a  height  of  200  km  the  mean  free  path  is  over 
200  m  and  free-molecule  flow  will  probably  be  experienced  by  all  satellites  so 
far  launched.  Even  at  175  km,  where  the  mean  free  path  is  over  100  m,  only  a 
few  satellites,  e.g.  the  large  carrier  rockets,  are  likely  to  have  entered  the 
transition  region,  where  collisions  between  incident  and  re-emitted  particles 
lead  to  a  reduction  in  the  drag  coefficient.  A  satellite  whose  greatest  dimen¬ 
sion  is  about  1  m  will  continue  in  free-molecule  flow  down  to  about  140  km. 

An  important  parameter  governing  the  drag  coefficient  of  a  body  in  free- 
molecule  flew  is  the  molecular  speed  ratio  s  defined  by 

speed  of  satellite  .  „ 

8  ...  - - - - - -  .  (2) 

most  probable  molecular  speed 

The  most  probable  molecular  speed  v  at  any  altitude  is  given  by 


where  T  is  the  atmospheric  temperature,  H  the  molecular  weight  and  R  the  gas 
constant.  vm  is  related  to  the  scale  height  H  by  the  equation 

v^  =  V&g  ,  (4) 

g  being  the  local  value  of  the  acceleration  due  to  gravity.  Using  values  of  H 

from*  the  COSPAR  International  Reference  Atmosphere  (19&4),  the  most  probable 

molecular  speed  is  found  to  increase  from  about  0.7  kn/sec  at  a  height  of  150  km 

to  around  1.5  kn/sec  at  fQO  km,  whereas  the  satellite’s  speed  is  always  greater 

than  7*5  kn/sec  at  these  altitudes.  Consequently  the  molecular  speed  ratio  at 

perigee  always  exceeds  5  for  perigee  heights  belcw  700  km  and  the  randan  thermal 

motion  of  the  atmospheric  molecules  can  be  neglected,  i.e.  the  flow  is  hyper- 
5 

thermal  •  Far  a  satellite  'with  perigee  at  1500  km,  however,  the  molecular  speed 
ratio  is  less  than  4  and  the  thermal  motion  cannot  be  ignored. 

3  LICIECUIAR  REFLEXION  AT  A  SCIJD  SURFACE 

3*1  The  accommodation  coefficient 

The  degree  of  energy  transfer  between  a  gas  molecule  and  a  surface  is 
usually  expressed  in  terms  of  the  accomodation  coefficient  a,  which  is  defined 
as  the  ratio  of  the  energy  change  experienced  by  the  Impinging  molecules  to  the 
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maximum  energy  change  that  could  take  place.  If  the  average  kinetic  energies  of 
the  incident  molecules  and  the  re-emitted  molecules  are  denoted  by  E^  and  Ey 
respectively,  and  the  average  kinetic  energy  of  a  molecule  re-emitted  with  ft 

velocity  corresponding  to  the  surface  temperature  is  E  ,  then 

s 


a 


(5) 


At  satellite  speeds  the  inoident  molecules  are  almost  uni-directional  and 
they  impinge  on  the  surface  with  energies  in  the  range  0.2  eV  (hydrogen  atom  at 
6  knysec)  to  12  eV  (oxygen  atom  at  8.5  ley 'sec) .  Consequently  we  are  interested 
in  incident  energies  which  are  high  compared  v/ith  the  energy  at  whicn  molecules 
are  likely  to  be  trapped  at  the  surface  by  physisorption,  but  still  v/ell  below 
the  level  at  which  appreciable  sputtering  will  occur.  So  far  no  experimental 
measurements  of  accommodation  coefficients  have  been  made  in  thi3  energy  range  . 
Over  the  past  few  years,  however,  a  nunber  of  theoretical  studies  have  been 
performed  and  these  liave  indicated  the  main  factors  affecting  the  energy 
accommodation. 


5.2  Theoretical  values  of  a 

The  theoretical  investigation  of  accommodation  coefficients  is  extremely 
difficult  because  cf  the  large  number  of  fundamental  assumptions  required  to 
make  progress.  These  assumptions  concern  the  motion  of  the  gas  atom  and  the 
configuration  and  motion  of  the  surface  atoms.  It  is  also  usual  to  assume 
that  the  surface  is  clean. 


In  terms  of  classical  theory  the  simplest  possible  model  is  obtained  by 
considering  the  incident  and  surface  atoms  to  be  smooth  and  perfectly  elastic 
spheres.  The  energy  exchange  i3  then  independent  of  gas  temperature  and  is 
greatest  in  a  head-on  collision  for  which  we  have  the  elementary  result 


a 


La 

(i  +  n)2  ’ 


(6) 


where  p  is  the  ratio  of  the  mast,  of  the  incident  gas  atom  to  the  mass  of  the 

surface  atom*  The  value  given  by  (6)  is  valid  for  p<1,  when  only  a  single 

collision  takes  place.  For  p  >  1  the  incident  gas  atom  is  still  travelling 

towards  the  surface  after  a  single  collision  and  ihe  accommodation  coefficient 

cannot  be  estimated  without  considering  the  motion  of  the  surface  atoms.  Far 

an  oblique  collision,  when  the  velocity  of  the  incident  atom  is  inclined  at 

2 

an  angle  <j>  to  the  line  of  centres  at  impact,  a  factor  cos  <j>  is  introduced  on 
the  right-liar  id  side  of  (6). 
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7-9 

Using  classical  mechanics  Goodman  lias  recently  made  an  extensive  study 
of  accommodation  coefficients  assuming  the  surface  to  be  a  simple  cubic  lattice. 
The  incictent  atom  is  assumed  to  impinge  directly  onto  one  cf  the  surface  atoms 
and  the  latter  are  considered  to  be  initially  at  rest.  In  his  latest  paper^  a 
Horse  potential  is  used  to  represent  the  interaction  between  the  gas  atom  and  the 
surface  atom.  For  a  single  atom  the  general  form  of  the  variation  of  the 
accommodation  coefficient  with  incident  energy  is  independent  of  the  gas- 
solid  system  and  is  as  shown  in  Fig- 2-,  At  lav  temperatures  (i.e.  lav  incident 
energies)  the  gas  atom  is  trapped  and  the  energy  accommodation  is  complete. 

Above  a  certain  temperature  the  gas  atom  lias  sufficient  energy  to  escape  and 
the  energy  accommodation  is  incomplete.  After  reaching  a  minimum  the  accommoda¬ 
tion  coefficient  increases  towards  the  value  given  by  the  elastic-spheres  model 
as  the  incident  energy  becomes  large.  Using  the  calculated  values  of 
accommodation  coefficient  for  single  atoms,  Goodman  proceeds  to  find  the  'true' 
accommodation  coefficient  on  the  assumption  that  the  incident  beam  has  a 
one-dimensional  Maxwellian  velocity  distribution.  If  the  parameters  of  the 
Horse  potential  are  suitably  chosen,  the  resulting  values  of  accommodation 
coefficient  show  adequate  quantitative  agreement  with  the  experimental  data 
of  Thomas  and  Schofield  (mentioned  in  Ref.  5)  over  a  limited  range  of  gas 
temperature. 

The  high  energy  limit  given  by  equation  (6)  will  only  apply  when  a  single 
collision  takes  place.  It  appears  for  the  elastic-spheres  model  and  also  far 
other  one-dimensional  lattice  models  that  a  second  collision  does  net  occur 
until  p  exceeds  1.  TTith  a  three-dimensional  lattice  model,  however,  Goodman  has 
foind^  that  there  is  always  a  second  collision  when  p  exceeds  O.84.  Consequently 
the  accommodation  coefficient  is  reduced  below  the  value  given  by  (6)  far 
0.84  <  p  <  1  and  it  becomes  small  as  p  tends  to  infinity.  Poor  accommodation 
for  large  values  of  p  has  also  been  noted  by  Gilbey10  for  a  continuum  model  cf 
the  surface.  Thi3  phenomenon  is  due  to  the  fact  that  a  tiireo-diaensional 
lattice  possesses  rigidity  and  always  completely  restores  itself  after  a  series 
of  collisions,  whereas  a  one-dimensional  lattice  docs  not. 

It  is  informative  to  consider  the  high-energy  limit  obtained  when  a 

quantum  theoretical  approach  is  used.  Again  an  incident  gas  atom  is  assumed 

to  impinge  directly  onto  a  lattice  atesw  It  was  originally  believed  that 

first-order  quantum-mechanical  perturbation  theories  were  inadequate  far  the 

10 

treatment  of  high  energy'  collisions.  Gilbey  has  recently  re-examined 
Devonshire’s  theory  and  shown  that  it  can  successfully  predict  multi-phonon 
transitions  and  should,  therefore,  predict  correct  values  of  accccmodation 
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coefficients  for  both  small  and  large  energy  transfers,  provided  the  accommoda¬ 
tion  coefficient  is  itself  small.  When  the  incident  energy  of  the  gas  atom  is 
sufficiently  high,  he  found  that  Devonshire's  result  reduces  to  a  s  4p.  Thi3 
is  an  approximation  to  (6)  valid  for  small  values  of  p« 

So  far  we  have  only  mentioned  theories  which  assume  the  incident  gas  at  cm 
to  impinge  directly  on  a  lattice  atom.  Although  this  assumption  is  useful  in 
initial  investigations  it  is  difficult  to  tell  how  the  predicted  values  will 
differ  frcm  the  values  which  occur  in  practice.  It  seeins  that  only  three 
attempts  have  been  made  to  consider  collisions  in  which  the  gas  atom  does  not 
approach  and  leave  the  solid,  along  a  normal  to  its  s'lrface.  In  Baule's  original 
treatment  the  surface  was  represented  by  an  oscillating  cubic  lattice. 
Collisions  were  based  on  the  elastic— spheres  model  and  the  energy  transfer  was 
determined  by  averaging  over  all  angles  of  incidence  <j>%  The  result,  obtained 
assuming  only  a  single  collision  with  the  surface,  gave 

2H 

“  =  T* - \Z  * 

(1  +  p) 

which  is  half  the  value  predicted  by  equation  (6). 

In  practice,  neither  equation  (6)  nor  equation  (7)  will  give  the  correct 

result,  and  the  true  value  of  a  will  probably  lie  somewhere  between  these  two 
1 2 

extremes.  $pik  has  also  attempted  to  calculate  average  values  of  accctmoda- 

tian  coefficients  using  the  elastic-spheres  model.  He  made  some  allowance 

for  the  fact  that  all  impacts  are  not  head-on  collisions  and  determined  an 

average  accommodation  coefficient.  His  computed  values  of  a,  which  are 

slightly  less  than  those  predicted  by  equation  (1),  are  shown  in  Pig.  3, 

together  with  these  from  other  theoretical  models,  hare  recently,  Oman,  Bogan, 
13 

Wei3er  and  Li  have  used  a  detailed  classical  model  in  which  gas  atom 

trajectories  are  not  confined  to  head-on  collisions.  The  results  of  sample 

14 

calculations  for  helium,  neon  and  argon  on  nickel  are  shown  in  Pig. 3  far 

o 

atomic  beams  incident  at  15  or  kb  to  the  surface  normal.  These  results 
suggest  that  a  value  between  the  two  curves  representing  equations  (6)  and  (7) 
is  of  the  right  order  at  energies  of  a  few  electron  volts. 

3  3  Distribution  of  reflected  molecules 

As  with  the  accommodation  coefficient,  no  measurements  of  the  distribution 
of  reflected  molecules  have  been  made  in  the  appropriate  energy  range.  Experi¬ 
mental  evidence  *hat  dees  exist  i3  to  sane  extent  conflicting,  and,  in  contrast 
to  the  accommodation  coefficient,  there  axe  no  theories  which  can  be  used  to 
assist  in  the  reconciliation  of  the  results. 
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The  extreme  mechanisms  of  re-emission  are  generally  regarded  as  specular 
and  diffuse.  Specular  reflexion  occurs  when  the  angle  of  reflexion  of  a  particle 
is  equal  ' o  the  angle  of  incidence.  Diffuse  re-emission  takes  place  according  to 
the  Knudsen  cosine  law,  which  states  that  the  number  of  molecules  emitted  between 
the  angles  0  and  0  +  6 0  from  the  normal  to  the  surface  is  proportional  to 
cos  0  60. 

At  thermal  energies  both  diffuse  re-emission  and  lobate  reflexion  near  the 

specular  ray  have  been  obtained  in  laboratory  experiments  with  molecular  beams. 

The  early  investigations  using  this  technique  were  principally  concerned  with 

scattering  by  cleaved  alkali  halides  and  this  work  has  recently  been  discussed 
15 

by  Hurlbut  .  For  these  surfaces  quasi-specular  reflexion  can  occur  at  all 

angles  of  incidence.  The  scattering  of  deuterium  and  helium  from  polycrystalline 

16 

platinum  has  been  studied  by  Data,  Koore  and  Taylor  ,  who  found  that  the  angalar 
distribution  conformed  to  two  patterns,  depending  on  the  temperature  of  the 
surface.  On  a  surface  at  roam  temperature  the  distributions  were  diffuse,  whereas 
on  a  high  temperature  surface  (about  1900°C)  they  were  remarkably  specular.  Far 
both  gases  the  beam  reflected  from  the  hot  surface  was  found  to  lie  mare  nearly 
along  the  specular  ray  as  the  beam  temperature  was  increased,  an  effect  that  has 
been  noted  by  several  investigators.  The  reflexion  of  low  energy  beams  of  rare 
gases,  nitrogen  and  oxygen  frcro  polycrystalline  platinum  and  nickel  has  been 
studied  recently  by  Hinchen  and  Foley  .  Again  diffuse  re-emission  was  observed 
when  the  surfaces  were  at  room  temperature.  Far  argon  on  platinum  a  change 
from  diffuse  re-emission  to  lobate  reflexion  was  found  a3  the  surface  was 
heated,  the  transition  occurring  in  the  region  of  500 °C.  The  lobe  maxima  were 
found  to  be  displaced  from  the  specular  ray  t wards  the  normal,  the  displacement 
being  greater  far  incident  angles  nearer  the  surface.  The  difference  in  the 
observed  distributions  between  the  hot  and  cold  surfaces  is  most  probably  due 
to  the  presence  of  adsorbed  gaseous  contaminants  at  low  temperatures. 

3*4  Application  to  satellites 

For  light  gases  on  clean  surfaces  both  theory  and  the  limited  experimental 
evidence  now  available  suggest  that  low  values  of  the  accomodation  coefficient 
are  appropriate.  The  success  of  Goodman's  theory  in  predicting  the  correct 
dependence  of  the  accomodation  coefficient  on  various  parameters  cannot  be 
ignored,  and,  for  a  clean  surface,  the  accomodation  coefficient  will  probably 
vary  with  incident  energy  in  the  manner  shown  qualitatively  in  Fig. 2.  The 
accommodation  coefficient  at  the  surfaoe  of  a  satellite  can  probably  be  cat- 
puted  on  the  basis  of  the  liard-spheres  model  so  that  a  is  proportional  to 


2 

|i/(l  +  n)  far  (i  <  1.  Since  the  factor  of  proportionality  is  likely  to  lie 
between  2  and  A,  the  accommodation  coefficient  will  probably  have  a  value 
between  the  continuous  and  dashed  curves  shewn  in  Pig»3« 

Before  an  acccmnodatiau  coefficient  can  be  chosen  for  a  particular 
satellite  at  a  specific  altitude  it  is  neoessary  to  decide  on  an  average  value 
for  p*  The  main  e  cmospheric  constituents  are  Ng»  0^,  0,  He,  A  and  H.  Atomic 
oxygen  predominates  over  most  of  the  altitude  range  of  chief  interest 
(150  to  800  km)  vhile  atomic  hydrogen  is  only  present  in  significant  quantities 
above  5C0  km.  Strictly  the  drag  due  to  each  constituent  should  be  computed 
separately,  but  because  of  other  uncertainties  it  is  adequate  to  take  the  mass 
of  the  incident  molecules  equal  to  the  ambient  molecular  weight  if  p  is  less 
than  one.  The  variation  of  molecular  weight  with  altitude,  as  given  by  the 
COSPAR  International  Reference  Atmosphere,  19 &*-,  shewn  in  Fig.4  far  extreme 

levels  of  solar  activity  and  local  times  corresponding  to  maximum  and  minimum 
values  of  the  diurnal  density  variation. 

When  the  incident  energy  of  the  impinging  molecules  is  sufficiently  large 
far  the ’hard-spheres'  model  to  be  valid,  the  mass  of  the  surface  atoms  used  in 
calculating  p  will  depend  cm  the  nature  of  the  outermost  layer  of  atoms  on 
the  surface,  and  not  on  the  actual  material  used  in  constructing  the  vehicle. 

For  only  a  few  satellites  have  the  surface  materials  been  publicized.  Typical 
materials  used  an  the  outside  of  satellites  are  gold,  stainless  steel,  aluminium, 
plastic,  quartz,  synthetic  sapphire  and  various  paints.  Quartz  or  synthetic 
sapphire  is  used  to  protect  solar  cellB  from  radiation  damage,  while  paints  are 
used  to  vary  the  surface  characteristics  and  hence  to  control  the  temperature 
of  the  vehicle. 

It  i3  extremely  difficult  to  define  the  exact  state  of  a  solid  surface 
owing  to  the  case  with  v/hich  gases  are  adsorbed.  The  type  of  adsorption 
depends  on  the  nature  of  the  attractive  force  holding  the  adsorbed  gas  to 
the  surface.  Physical  adsorption  depends  on  van  der  V/aals '  forces  and  the  heat 
of  adsorption  is  relatively  low,  i.e.  a  few  kcal/mole.  Chemisorption  displays  a 
large  heat  of  adsorption,  implying  the  existence  of  a  chemical  bond.  In  contrast 
to  physisorption,  chemisorption  is  not  readily  reversible,  high  temperatures  and 
low  pressures  being  required  to  produce  desorption. 

Chemisorption  tokos  place  under  conditions  in  which  incipient  chemical 
reaction  between  gas  and  surface  might  be  expected.  For  example,  most  metals 
are  capable  of  reacting  with  oxygen  to  farm  oxides,  so  that  it  is  not  sur¬ 
prising  tiiat  oxygen  is  chemisorbed  by  all  metals  except  gold.  Although  an 
energy  of  activation  is  sometimes  necessary  before  chemisorption  occurs,  it  can 
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occur  extremely  rapidly,  especially  with  oxygen,  and  a  complete  monolayer  can 
be  f conned  on  a  metal  surface  in  a  few  seconds,  even  at  a  pressure  of  10~^  nm 
of  mercury  (which  corresponds  to  a  height  of  about  200  km). 

Since  oxygen  is  fairly  abundant  at  ground  level  and  atomic  oxygen  is 
usually  the  main  atmospheric  constituent  at  the  orbital  heights  of  interest, 
we  would  expect,  on  the  basis  of  the  foregoing  discussion,  that  the  surfaces 
of  all  metals  except  gold  would  be  covered  by  a  monolayer  of  adsorbed  oxygen, 
ii any  metals  would  possess  mare  than  a  monolayer  since  oxide  films  con  be  formed 
by  a  diffusion  mechanism.  Aluminium  is  known  to  farm  an  oxide  film  in  this 
manner,  while  stainless  steel  ewes  its  'stainless1  property  to  the  formation 
of  a  film  of  chromic  oxide.  It  is  reasonable  to  suppose  that  in  the  final  atomic 
arrangement  the  metallic  surfaoes  of  a  satellite  would  most  probably  present 
oxygen  atoms  to  the  ambient  atmospheric  particles;  only  gold-plated  surfaces 
are  likely  to  remain  uncontaminated.  Although  an  atom  will  occasionally  be 
desorbed,  adsorption  sites  are  unlikely  to  remain  vacant  far  long  in  the 
presence  of  energetic  atoms  of  oxygen. 

For  non-metals  the  nature  of  the  outermost  layer  of  atoms  is  not  so  easily 
defined.  Physisarbed  atoms  will  be  readily  desorbed  by  the  incident  atoms 
while  the  chemisorption  of  energetic  atoms  on  non-metals  is  virtually  unexplored. 
Sapphire  is  and,  since  the  ionic  radius  of  aluminium  is  small  compared 

with  that  of  oxygen,  a  gas  atom  is  unlikely  to  hit  aluminium.  With  quartz 
(Si  02) ,  a  free  space  near  a  silicon  atom  would  constitute  a  site  for  chemisorp¬ 
tion  of  an  oxygen  atom  and  would  be  unlikely  to  remain  unfilled  far  long.  Since 
quartz  possesses  a  somewhat  open  structure,  an  occasional  hit  on  a  silicon  atom 
cannot  be  ruled  out.  The  atomic  weights  of  aluminium  and  silicon  are,  however, 
not  so  very  different  from  oxygen,  and  the  small  uncertainty  in  (i  will  not 
significantly  contribute  to  the  already  large  uncertainty  in  the  accommodation 
coefficient. 

The  outermost  atoms  of  a  painted  surface  are  even  mare  difficult  to  specify. 
Paints  consist  of  inorganic  pigments  and  binders  of  silicone^  acrylic  or  epoxy 
resins.  These  resins  are  synthetic  polymers  with  organic  groups  on  the  polymer 
chain.  Since  the  organic  groups  contain  hydrogen,  the  surface  atoms  may  some¬ 
times  be  lighter  than  the  incident  atoms,  which  would  then  undergo  a  second 
collision.  At  present  it  i3  impossible  to  specify  even  an  approximate  value  of 
a  under  these  circumstances.  Fortunately  it  is  not  usually  necessary  to  paint 
mare  than  20  to  3 Q*  of  a  satellite's  surface  far  thermal  control  purposes. 
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A  topic  which  should  be  mentioned  at  this  stage  is  sputtering.  A  few 

years  ago  it  was  common  practice  to  discuss  threshold  energies  below  which 

sputtering  did  not  occur,  values  of  50-60  eV  being  quoted  far  most  materials. 

With  the  mare  accurate  measuring  techniques  now  available,  very  small  yields 

18 

have  been  measured  at  lower  energies.  Extrapolation  of  these  results  suggests 

that  sputtering  is  insignificant  at  energies  below  20-25  eV.  It  is  therefore 

reasonable  to  conclude  that  sputtering  will  have  a  negligible  effect  an  drag, 

both  directly  through  the  removal  of  momentum  by  sputtered  material,  and 

indirectly  through  the  modification  of  p  by  the  removal  of  adsorbed  atoms. 

19 

McKeown  and  his  collaborators  claim  to  have  measured  sputtering  an  two 

Discoverer  satellites  at  a  height  of  about  200  km.  The  erosion  rates  were 

found  to  be  0.1  ±  0.05  X/day  far  gold  and  0.15  ±  0.05  2/aay  for  silver,  corres- 

-6  -6 

ponding  to  sputtering  yields  of  1C~  atoms/molecule  and  2  x  1C  atcms/molecule 

20 

respectively.  These  values  have  recently  been  disputed  ,  however. 


4  DRAG  OF  SIliPlE  SHARES  IN'  HYIERTHERIIAL  EREE-MOIECUIE  FLOW 


From  equation  (5)  the  ratio  of  the  speed  of  a  re-emitted  molecule  vr  to 
Uie  speed  of  an  incident  molecule  v^  is  related  to  the  thermal  accommodation 
coefficient  by  the  equations 


and 


J 


(8) 


where  T^  is  the  kinetic  temperature  of  incidence  and  Tw  is  the  surface  tempera¬ 
ture.  Since  we  are  unable  to  make  any  statement  about  the  variation  of  a  with 
the  angle  of  incidence  of  a  surface  element,  all  reflected  molecules  are 
assumed  to  have  the  same  value  of  vr»  After  making  tills  assumption,  there  is 
little  point  in  taking  into  account  the  random  thermal  motion  for  molecular 
speed  ratios  above  5* 

The  drag  of  satellites  in  hyperthermal  free-molecule  flow  has  been  studied 
by  Schamberg^  and  his  results  far  certain  simple  shapes  are  quoted  in  Table  1. 
The  drag  coefficient  is  defined  by 


°D  * 


D 


ipV2S 


(9) 
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where  D  is  the  drag  acting  in  the  opposite  direction  to  the  satellite's  velocity 
vector  V,  p  is  the  air  density  and  5  is  a  representative  cross-sectional  area. 
The  ratio  v^/v^  denoted  by  r.  Since  T is  very  small,  r  is  approximately 
equal  to  VI  -  a* 

Table  1 

DRAG  IN  HYEERTHERL1AL  FREE-MOIECUIE  FLOW 


Shape 


Plat  plate 

(normal  to  flow) 

Plat  plate  at 
incidence  6 

Sphere 

Cylinder  perpendicular 
to  flow 

Cone  of  semi  vertex 
angle  i|r  with  vertex 
forwards  and  axis 
parallel  to  flow 


Drag  coefficient 

(based  on  projected  area  perpendicular  to  direction 
of  motion) 

Diffuse  re-emis3ion 

Accommodated  specular  reflexion 

2(1  +  j  r) 

2(1  +  r) 

2(1  +  -j  r  sin  6) 

2(1  -  r  cos  2  6) 

2(1  ♦  |  r) 

2 

2(1  +  f  r) 

2(1  *ir) 

2(1  +  |  r  sin  V) 

X— N 

CM 

01 

O 

O 

u 

1 

w— 

CM 

The  variation  of  the  drag  coefficient  with  a  is  shown  in  Pigs.  5  to  7  far 
a  flat  plate,  a  sphere  and  a  circular  cylinder  of  infinite  length  respectively. 

Cp  was  evaluated  assuming  T^/T^  =  0.006.  Pig.5»  which  gives  the  drag  coefficient 
of  a  flat  plate  based  on  the  projected  area  perpendicular  to  the  direction  of 
motion  for  various  angles  of  incidence,  serves  as  a  reminder  that  Cp  would  not 
necessarily  exceed  2  if  there  were  appreciable  reflexion  of  molecules  along  the 
specular  ray.  A  flat  plate  is  not  a  practical  shape,  however,  and  Pigs. 6  and  7 
give  a  fairer  impression  of  the  likely  value  of  Cp. 


5  CRAG  OF  A  SPINNING  CIRCULAR  CYLINDER  IN  HYTSRTHERliAL  PREE-MOLBCULE  PLOW 


Most  satellites  used  for  density  deternunation  are  unstabilized,  so  that 
their  attitudes  are  variable.  It  is  known  that  a  spinning  satellite  will 

eventually  rotate  about  its  axis  of  maximum  moment  of  inertia  if  it  is  either 

21  22 
acted  upon  by  small  external  torques  or  experiences  energy  dissipation  internally  . 

Hence,  to  evaluate  drag,  a  ncn-spherical  satellite  is  normally  assumed  to  rotate 

abort  its  axis  of  maximum  moment  of  inertia  after  it  has  been  in  orbit  for  a  few 
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days.  Unfortunately  the  position  of  the  axis  of  rotation  is,  in  general, 
unknown.  The  extreme  modes  of  rotation  for  a  long  thin  satellite  spinning 
about  its  axis  of  maximum  moment  of  inertia  are  (a)  spinning  propellerwise 
and  (b)  tumbling  end-over-end.  Since  the  axis  of  rotation  may  be  inclined 
at  any  angle  to  the  direction  of  motion,  it  is  best  to  take  for  the  drag  the 
mean  of  its  values  under  modes  of  rotation  (a)  and  (b). 

Many  spinning  satellites  are  approximately  circular  cylinders  and  this 
is  the  only  shape  which  need  be  considered  in  detail.  The  drag  wider  mode  (a) 
has  already  been  given  in  the  previous  section;  it  remains  to  oonsider  the 
drag  under  mode  (b). 

5 

Using  the  equations  given  by  Schamberg  far  the  farces  acting  on  an 
element  of  area,  the  drag  coefficient  of  a  circular  cylinder  tumbling  end-over- 
end  is  readily  shown  to  be  given  by 


°D  = 


(10) 


for  diffuse  re-emission.  Here  C  is  the  length  of  the  cylinder  and  d  is  its 
diameter.  The  reference  area,  which  i3  taken  as  the  mean  of  the  projected 
area  perpendicular  to  the  direction  of  motion,  iB  given  by 

3  -  *  x)  •  (11) 


The  drag  coefficient  given  by  equation  (10)  is  almost  independent  of  the 
length/diameter  ratio  as  can  be  seen  from  Fig.8.  It  is  worth  noting  that  the 
drag  coefficient  of  a  rotating  body  is  close  to  the  value  far  a  sphere  if  the 
mean  projected  area  is  taken  as  the  reference  area. 

6  DRAG  AT  LG?  MOLECULAR  SHSED  RATIOS 

As  the  altitude  increases  above  about  800  km  the  most  probable  molecular 
speed  bee  ernes  so  large  that  the  random  thermal  motion  of  the  atmospheric 
molecules  must  be  considered.  The  drag  produced  by  impinging  molecules  is 
readily  calculated,  but,  unless  specular  reflexion  occurs,  the  drag  due  to  the 
re-emitted  molecules  is  not  so  readily  obtained. 

Here  we  consider  only  spheres,  for  which  the  drag  due  to  the  incident 
molecules  only  is  the  same  as  the  total  drag  with  specular  reflexion.  If  i 
refers  to  incident  molecules  and  S  to  specular  reflexion,  the  corresponding 
drag  coefficients  are  given  by 


14 


oc 


2 

\  ‘  °D,  •  \1+7-*k)  erf(s>  *  axp  (-"2)  •  <12> 

1  23  2L 

When  considering  drag  with  diffuse  re-emission  most  authors  '  assume  that 
the  most  probable  re-emission  speed  vr  is  constant  over  the  entire  surfaoe. 

The  contribution  of  the  re-emitted  moleoules  to  the  drag  coefficient  is  then 
given  by  2/V3sr  so  that  the  total  drag  is 

°D  '  •  <15) 

Although  the  assumption  of  constant  re-emission  speed  is  valid  far  complete 
accomodation,  it  does  not  apply  for  poor  accomodation;  so,  for  molecular 
speed  ratios  above  one,  the  molecular  speed  ratio  of  re-emission  will  be  much 
larger  at  the  rear  than  at  the  front* 

With  present  knowledge  of  gas-Burfaoe  interactions,  detailed  calculations 

allowing  far  variations  in  the  re-emission  speed  over  the  surface  are  not 

justified.  Far  accommodation  coefficients  near  zero  and  molecular  speed  ratios 

greater  than  one,  it  is  at  present  preferable  to  assume  one  re-emission  speed 

ratio  s  far  all  surfaoe  elements  at  the  front  and  another  re-emission  speed 
r 

ratio  s'  far  all  surface  elements  at  the  rear.  It  can  then  be  shown  that  the 
r 

drag  coefficient  due  to  diffuse  re-emission  is  given  by 


so  that  the  total  drag  coefficient  is 

°D  *  \*°B.  •  ‘15> 

x  CL 

We  can  calculate  what  is  probably  the  largest  possible  value  of  by  taking 

s  a  1  and  8 '  s  s.  This  means  that  we  are  assuming  the  moleoules  incident 
r  r 

on  the  front  are  re-emitted  with  a  speed  equal  to  the  satellite's  speed,  and 
the  moleoules  incident  on  the  rear  are  re-emitted  at  a  speed  corresponding 
to  the  asfcient  atmoepherio  temperature.  Fig. 9  indicates  that  a  practical  upper 
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limit  to  the  drag  coefficient  ie  probably  boat  estimated  by  drawing  a  curve 
from  the  hyperthermal  limit  at  s  s  6  to  the  upper  limit  for  the  drag  at  s  =  1. 

The  lowest  curve  in  Fig. 9  gives  the  drag  coefficient  of  a  sphere  for 
either  specular  reflexion  or  diffuse  re-omiesion  with  T ^/T  =  0,  Tr  being  the 
temperature  of  re-emission.  The  drag  with  diffuse  re-emission  and  complete 
accommodation  can  be  obtained  for  any  value  of  J T ^/T  by  linear  interpolation 
between  the  three  lowest  curves.  The  drag  coefficient  for  diffuse  re-emission 
and  incomplete  accommodation  can  be  obtained  by  interpolating  between  the 
praotical  upper  limit  discussed  above  and  the  appropriate  value  for  complete 
accommodation,  bearing  in  mind  the  variation  of  with  a  far  hyperthermal 
free-molecule  flew  given  in  Fig. 6. 

For  values  of  s  near  and  belcw  unity  the  most  probable  molecular  speed 
is  close  to  the  local  escape  speed  and  the  assumption  of  a  Maxwellian 
velocity  distribution  in  the  free  stream  is  no  longer  Justified.  Far  a 
circular  orbit  the  drag  coefficient  should  not  be  taken  frem  Fig. 9  far  mole¬ 
cular  speed  ratios  less  than  about  1*5. 

7  CHARGE  DRAG 

A  satellite  moving  in  the  upper  atmosphere  can  acquire  an  electric  charge 

and  there  will  then  be  a  contribution  to  the  drag  from  the  interaction  of  the 

vehicle  with  the  ambient  ions.  Several  analyses  of  charge  drag  have  been  mads 

25  26 

and  many  conflicting  results  obtained.  On  the  basis  of  the  present  evidence  , 
charge  drag  appears  to  be  negligible  up  to  at  least  the  altitude  where  helium 
becomes  the  main  constituent.  This  occurs  near  550  km  at  night  when  solar 
activity  is  low  and  at  much  higher  altitudes  at  other  times.  Charge  drag  may 
still  be  negligible  at  much  greater  altitudes,  however. 

27 

An  experimental  investigation  of  the  drag  of  six  tin  dipoles,  34  cm  long 
and  0.043  cm  in  diameter,  at  an  altitude  of  about  3000  km  has  shown  that  the 
effect  of  charge  drag  on  small  satellites  at  this  altitude  is  not  catastrophic, 
as  has  sometimes  been  predicted,  although  it  may  exceed  the  neutral  particle 
drag.  The  smallness  of  the  electric  drag  1ms  been  confirmed  by  observation  of 
the  West  Fold  test  belt23. 

8  DISCUSSION 

In  the  light  of  the  discussion  cm  surface  interaction  in  section  3#  it  is 
evident  that  there  is  still  considerable  uncertainty  in  the  estimation  of 
satellite  drag  coefficients.  There  is  new  ample  evidence,  both  theoretical 
and  experimental,  that  lew  values  of  accommodation  coefficient  are  appropriate 


for  light  gases  incident  on  heavy  surface  atoms  and  a  can  be  taken  to  increase 
fran  zero  for  a  mass  ratio  p  of  zero  to  a  value  near  unity  for  a  mass  ratio  of 
1.  We  have  seen  that  the  mass  ratio  depends  on  the  outermost  layer  of  atoms 
on  the  surface  and  on  the  ambient  molecular  weight.  Most  satellite  surfaces, 
unless  they  are  either  gold  plated  or  painted,  will  present  mainly  atoms  of 
oxygen  to  the  incident  gas  so  that  p  will  be  less  than  1  when  the  ambient 
molecular  weight  is  less  tlian  about  16.  In  the  important  region  where  the 
molecular  weight  exceeds  16  it  is  probably  best,  at  present,  to  take  the  value 
of  a  appropriate  tc  p  =  1  for  the  majority  of  surfaces.  This  is  because  high 
energy  collisions  take  place  practically  atom  with  atom.  Also,  since  the 
dissociation  energies  of  Og  and  Ng  molecules  are  5*08  and  9.76  eV  respectively, 
there  will  usually  be  seme  dissociation  at  low  altitudes  and  this  will  have 
the  effect  of  increasing  the  accommodation  coefficient. 

The  exact  distribution  of  the  reflected  molecules  is  unknown  and  it  is 
probably  best  to  assume  diffuse  re-emi3sicm  until  further  experimental  evidence 
is  available.  Fortunately,  for  accommodation  coefficients  near  unity,  the  drag 
coefficient  is  not  particularly  sensitive  to  the  distribution  of  the  reflected 
molecules  and  little  error  will  ensue  from  assuming  an  incorrect  distribution. 

As  the  accommodation  ooeffioient  decreases,  however,  the  drag  becomes  mare 
sensitive  to  the  distribution  of  the  reflected  molecules. 

In  the  absence  of  a  proper  understanding  of  gas-surface  interactions,  it 
is  only  appropriate  to  use  an  elementary  model  to  calculate  drag.  Far  satellites 
below  700  km  the  molecular  speed  ratio  exceeds  5  and  the  random  thermal  motion 
of  the  atmospheric  molecules  can  be  ignored.  The  drag  of  various  simple  shapes 
can  be  obtained  from  the  formulae  in  Table  1  or  from  Figs.  5  to  8.  For  purposes 
of  illustration  the  variation  of  the  drag  coefficient  of  a  sphere  with  altitude 
is  shewn  in  Fig.  10  for  three  variations  of  a  with  p.  The  four  curves  for  each 
variation  of  a  correspond  to  the  values  of  molecular  weight  shown  in  Fig. 4. 

The  accommodation  coefficient  at  the  surface  of  most  satellites  probably 
exceeds  0.8  at  heights  up  to  about  400  km,  at  all  times  of  day  and  far  all  levels 
of  solar  activity.  For  a  sphere  the  smallest  value  of  the  drag  coefficient  is 
then  2.07  and  occurs  when  a  is  unity.  A  drag  coefficient  of  2.2,  which  is  the 
value  normally  used  in  density  determinations,  corresponds  to  a  &  0.95*  it 
must  be  emphasized  that  this  value  may  well  be  too  low,  perhaps  by  as  much  as 
10  per  cent,  especially  if  a  decreases  appreciably  as  the  angle  between  the 
incident  beam  and  the  surface  normal  increases.  Nevertheless  the  possible 
range  of  values  of  CL.  givon  by  the  theory  discussed  above  is  insufficient  to 
explain  7  the  large  discrepancy  (a  faotar  of  2)  between  the  density  values 
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measured"1^  by  ionisation  gauges  aboexd  Explorer  17  (1 963-0 9A)  at  a  height  of 
280  km  and  those  obtained  from  orbital  studies. 

At  altitudes  above  400  km  the  mass  ratio  (i  can  be  quite  small,  particu¬ 
larly  a+  times  of  low  solar  activity,  and  a  somewhat  higher  value  of  is  then 
appropriate.  The  drag  coefficient  of  any  particular  satellite  can  only  be 
obtained  after  making  a  judicious  choice  of  the  various  parameters  involved. 

The  recent  discussion  of  satellite  drag  by  Izakov^  is  oversimplified,  since 
the  degree  of  energy  transfer  is  assumed  to  be  the  same  at  all  heights. 

For  satellites  above  700  to  800  km  it  is  neoessary  to  take  the  thermal 
motion  into  account  and  the  remarks  made  in  section  6  must  be  borne  in  mind. 

The  drag  coefficients  of  spheres  in  high-altitude  oiroular  orbits  are  given  in 
Table  2  for  various  values  of  exospheric  temperature  and  two  values  of  molecular 
weight.  The  values  of  were  obtained  from  Fig.  9  by  interpolation,  assuming 
a  =  3«6p/(l  +  |i)  Cq  does  not  vary  appreciably  if  the  numerical  factor  in 
this  equation  is  changed  to  3. 2  or  4»  A  temperature  of  JCo\  corresponds  to 

Table  2 

DRAG  COEFFICIENTS  OF  SPHERES  IN  IIIGH-AMITUDE  CIRCULAR  ORBITS 


Molecular  weight 

1 

4 

Temperature  (°K) 

700 

1000 

1000 

1500 

2000 

Height  (kn) 

■ 

1C00 

3.4 

3.6 

2.8 

2000 

3.5 

3.7 

2.8 

2.9 

ESI 

3000 

3.5 

3.8 

3.0 

51 

32 

very  lew  solar  activity  ,  when  atomic  hydrogen  is  likely  to  be  the  dominant 
constituent  above  about  800  km.  A  temperature  of  2000 corresponds  to  very 
high  solar  aotivity  when  tliere  is  a  strong  possibility  that  helium  will  be  the 
dominant  constituent  above  about  1300  km.  The  drag  coefficients  given  in 
Table  2  apply  if  the  surface  a tans  have  a  molecular  weight  of  16.  When  hydrogen 
is  a  major  constituent  these  values  of  serve  as  an  upper  limit,  since  far 
seme  surfaces  the  outmost  stems  may  not  then  be  oxygen. 

9  CONCLUSIONS 

Between  heights  of  140  and  400  ion  at  times  of  low  solar  activity,  or 
heights  of  1i*0  and  600  km  at  times  of  high  solar  aotivity,  the  drag  coefficient 
is  almost  independent  of  height,  and  at  present  for  most  satellites,  because  of 
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the  laok  of  deoisive  evidence,  there  is  perhaps  not  sufficient  reason  to  abandon 
the  value  of  2.2  whioh  has  been  widely  used  in  reoent  years  for  the  drag  co¬ 
efficient.  It  must  be  reoognized,  however,  that  this  value  is  subject  to  some 
uncertainty  and  may  be  too  lew,  perhaps  by  as  much  as  10  per  cent. 

At  heights  above  400  km  (low  solar  activity)  or  600  km  (high  solar 
activity)  the  drag  coefficient  increases  with  height,  since  both  the  energy 
transfer  and  the  molecular  speed  ratio  decrease  as  the  molecular  weight  of  the 
atmosphere  decreases.  It  should  be  noted  that  the  uncertainty  in  the  accommoda¬ 
tion  coefficient  is  less  when  the  molecular  weight  of  the  atmosphere  is  low. 

This  leads  to  only  a  small  spread  in  the  possible  values  of  drag  coefficient 
in  the  700-800  km  region  at  times  of  low  solar  aotivity,  as  can  be  seen  from 
Pig.10;  under  these  conditions  the  drag  coefficient  can  be  predicted  to  within 
a  few  per  cent. 
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SYMBOIS  (Cental 

X  mean  free  path  near  front  surface 

X  mean  free  path  in  free  stream 
o 

p  ratio  of  mass  of  gas  atom  to  mass  of  surface  atom 

p  air  density 

l  see  seoticn  3*2 

t  semi  vertex  angle  of  cone 
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FIG.  I  MEAN  FREE  PATH  FOR  HEIGHTS  BETWEEN  IOO  AND 
250  km,  AS  GIVEN  BY  THE  U.S.  STANDARD  ATMOSPHERE  ( 1962  ) 
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FIG.  2  FORM  OF  THE  VARIATION  OF  ACCOMMODATION 
COEFFICIENT  WITH  INCIDENT  ENERGY (aFTER  GOODMAN9) 
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FIG.  3  VARIATION  OF  ACCOMMODATION  COEFFICIENT 
WITH  MASS  RATIO  /A,  AS  GIVEN  BY  VARIOUS  THEORIES 
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FIG.  4  EXTREME  VALUES  OF  THE  MOLECULAR  WEIGHT  OF 
THE  ATMOSPHERE  AS  GIVEN  BY  CIRA  (l964) 


FIG. 5  DRAG  COEFFICIENT  OF  A  FLAT  PLATE  IN  HYPERTHERMAL 
FREE- MOLECULE  FLOW.  THE  DRAG  COEFFICIENT  IS  BASED  ON 
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FIG.  6  DRAG  COEFFICIENT  OF  A  SPHERE  IN  HYPERTHERMAL 
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FIG. 7  DRAG  COEFFICIENT  OF  AN  INFINITE  CIRCULAR  CYLINDER 
WITH  AXIS  PERPENDICULAR  TO  DIRECTION  OF  MOTION,  IN 
HYPERTHERMAL  FREE  -  MOLECULE  FLOW.(I*  -O  •  006) 
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FIG. 8  DRAG  COEFFICIENT  OF  A  CIRCULAR  CYLINDER  TUMBLING 
END-OVER -END,  WITH  DIFFUSE  RE- EMISSION.  THE  DRAG 
COEFFICIENT  IS  BASED  ON  THE  AVERAGE  PROJECTED 
AREA  PERPENDICULAR  TO  THE  DIRECTION  OF  MOTION. 
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FIG.  9  VARIATION  OF  THE  DRAG  COEFFICIENT  OF  A  SPHERE 


WITH  THE  MOLECULAR  SPEED  RATIO 


Fig. 10 


SPA/P/ 14  17 


200  400  6 00  800 

ALTITUDE  (ktn) 

FIG.  IO  VARIATION  OF  THE  DRAG  COEFFICIENT  OF  A 
SPHERE  WITH  ALTITUDE  ASSUMING  HYPERTHERMAL 

FREE  -  MOLECULE  FLOW 


